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Abstract. The coherent rise and fall of the quasar population is dis- 
cussed in terms of gas accretion onto massive black holes, governed by 
^ ■ the hierarchically growing environment. The rise is related to plentiful 

, accretion during the assemblage of the host galaxies; the fall to intermit- 

pLn ' tent accretion when these interact with companions in a group. The LFs 

^ , are computed out to z = 6, and are related to the mass distribution of 

^Sj ' relict BH found in local galaxies. The histories of the QS and of the star 

light are compared. 
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I The bright quasars exhibit remarkable permanence and remarkable variations: 

. the spectra (including the optical emission lines) are basically similar for indi- 

\ vidual QSs shining at redshifts in the full range z ~ — 5, throughout most 

Q> ■ of the universe life; meanwhile, the population undergoes substantial changes. 

\ These look well organized, as shown by surveys in the radio, the optical and 

Q^' the X-ray band (see Shaver et al. 1996, Osmer in this Volume, and references 

L \ therein). In fact, the bright QS population rises and falls steeply, culminating 

\^ • at2;~3ib0.5; fig. 1 represents this course in terms of optical light density, and 

\ compares it with the run of the star light density as given by Madau 1997. 

■ Over the first few Gyrs of the universe life the number of bright QSs grows, 

^ , and the luminosity functions mainly rise in normalization, a so-called negative 

density evolution. Instead, after z ~ 3 the optical LFs fall to nearly blend in 
with the local Seyfert 1 nuclei. Luminosity evolution apparently prevails, but 
' some positive density evolution also occurs; in addition, the optical LFs flatten 

toward us to comprise more bright objects than previously recognized (La Franca 
& Cristiani 1997; Goldschmidt k Miller 1997; Kohler et al. 1997). 

On the other hand, aimed images from HST and from the ground in optimal 
seeing conditions (Hutchings & Neff 1992, Disney et al. 1995, Hasinger et al. 
1996, Bahcall et al. 1997) extend to quasars out to 2 ~ 0.3 the evidence of a 
complex environment nailed down by Rafanelli et al. 1995 for the local Seyferts. 
In fact, up to 1/2 of the hosts mapped are found to be either engaged in galaxy 
interactions and in merging events, or to have close neighbors, even submerged 
within the host body. At a somewhat wider range, the QS environments are 
found to comprise some 10 — 20 galaxies on average (Fisher et al. 1996), the 
membership of a group. 

Thus the evidence points toward one basic engine, but with working regimes 
related to the surrounding structure. We shall discuss how such connection gives 
rise to non-monotonic evolution. 
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Figure 1. The history of the blue luminosity density produced by the 
QSOs with Lb > 10^^ erg is represented by the solid line (scale on 
the right). The points outline the UV light history of the stars after 
Madau 1997 (left scale). For the radio sources, see Dunlop 1997. 

2. Paradigms 

QSs as individual sources are held to be powered by gas accretion onto massive 
black holes. This had been widely argued to explain outputs up to L ~ 10^^ 
erg s~^ with large L/M and high apparent compactness, considering also the 
terminal stability of the BHs. The paradigm envisages extracting the power L 
from a mass accretion M, at distances r down to a few times the Schwarzschild 
radius rs and with radiative efficiencies up to r] = L/(?M ^ rg/r ^ 0.1. 

Direct evidence confirming the paradigm is now mounting (see Rees 1997), 
while the need for high rj is weakening. This is because relict BHs - expected to 
reside in normal galaxies keeping the record of the mass accreted during their 
whole career - are now found with number x masses exceeding the expectations 
under top efficiency (Kormendy k, Richstone 1995, Magorrian et al. 1997). 

But the paradigm has little to say - directly - concerning the rise and fall 
of the QS population. In fact, accretion drawing from an unlimited mass supply 
but self-limited by the BH radiation pressure in the Eddington regime yields 
luminosities that scale like Le — 10^^ M^jz/IO^Mq erg s~-^, and exponentiate 
on the short time scale -qtE — ^ 10^ yr, or last even less according to Haiman & 
Loeb 1997. Thus the activities of the individual QSs constitute just short flashes 
compared with the evolutionary times of a few Gyrs or with the population 
lifespan of ~ 12 Gyr. A complementary, coordinating agency is called for, to 
organize all those flashes into a coherent rise and fall over an Hubble time. 

Here enters the other paradigm, the hierarchical growth of structures which 
builds up the BH environment and the gas reservoir, and regulates the accretion. 
The gas may be held at bay on scales ^ 10^ pc in an axisymmetric gravitational 
potential where the angular momentum j is conserved. But such symmetry 
is broken during the host galaxy build Tip, when strongly asymmetric events 
of merging occur; these allow plentiful mass inflow. Subsequently, the hosts 
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stabilize but are enclosed in groups, where they interact with companion galaxies. 
The potential is again distorted giving rise to episodes of mass inflow, recurring 
but gradually petering out as groups are reshuffled into clusters. 

We maintain, following up CPV97, that in both dynamical regimes j at 
times is not conserved, providing a condition necessary for growing new BHs, or 
for refueling the old ones. The transitional mass from a large galaxy to a small 
group is around 5 10^^ Mq; in the hierarchical cosmogonies this corresponds to 
z ~ 3 lb 0.5, depending on cosmogonical and cosmological details. In any case, 
such values interestingly fall in the range where the bright QSs peak. 

3. Bimodal accretion 

Hierarchical cosmogony (see Peebles 1993) envisages larger and larger structures 
condensing out of gravitationally unstable density perturbations dominated by 
dark matter. In the critical universe the typical dark halos condensing and 
virializing at the epoch t, so attaining density contrasts p/pu{z) ^ 2 10^, scale up 
in mass after Mc oc for cold DM perturbations n slowly increases with 

M in the range ~ —2.5-. — 1.5. The halo growth may be visualized as a sequence 
of merging events of unequal, sometimes comparable, blocks. Thus if a typical 
rich cluster forms now, a small group formed at ^; ~ 2.5, and most galactic bulges 
formed before. But the actual condensed masses arc widely dispersed around 
Mc, correspondigly, their number rises prior to tc oc M*^"+"^)/^, but declines only 
slowly thereafter. Even in the adverse open cosmologies with Qo ^ 1 a similar 
trend is retained until the perturbations freeze out at 1 + z ~ l/^o- 

A subgalactic building block of DM with M ~ 10^° Mq allows a BH of 
nearly 10^ Mq to form, involving a baryon fraction around e ~ 10~^. The 
galaxy assemblage goes on by repeated, chaotic merging; the baryons lose an- 
gular momentum to the DM at a rate set by the number density of substruc- 
tures surviving in the gas and in the DM, and so approximately proportional 
to Puiz)- Inflow is plentiful and accretion only self-limited, to yield luminosities 
L ^ Le on Mbh oc e M where e may still vary with M and z as discussed below. 

In groups, many simulations (see Governato, Tozzi & Cavaliere 1996) have 
shown galaxy interactions to be frequent and strong. This is due to the high 
density of galaxies ng oc Pu{z) and to the low velocity dispersion V oc m'^ n)/i2 ^ 
still so close to the galaxian Vg as to allow dynamical resonance (such conditions 
no longer hold in clusters). Nearly grazing encounters occur frequently on the 
time scale ~ l/^'^g n^g and produce outright galaxy aggregations (Cava- 
liere & Menci 1997), cannibalism, or interactions strong enough to perturb the 
potential and cause j ^ const again. Simulations of single, aimed interactions 
(see Barnes & Hernquist 1991) show in detail how a sizeable fraction of the gas 
in both partners loses j and is driven toward the main galactic nucleus, down 
to a distance limited to now only by the computational dynamic range. 

So the BH bimodal fueling during host formation and interactions is unifled 
and coordinated by the hierarchical evolution of structures. This provides the 
necessary external condition for mass inflow, while ultimate acceptance by the 
BH is set by the radiation pressure. All that does not end the story, however. 

Initially, with copious inflow the stability of any intermediate structure 
against energy deposition conceivably limits the BH to Mbh oc (1 -|- z)"^'^ M^/^ 
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(HNR97), corresponding to e{M,z) oc (1 + z)^-^M^/^. At last, when external 
conditions allow only inflows M ^ 10^'^ Le/c^, one expects energy advection 
from the accretion disk down the BH horizon before the electrons share and 
radiate the ion energy. If so (but see Bisnovatyi-Kogan & Lovelace 1997), the 
residual emission should peak in radio and in hard X-rays (Di Matteo & Fabian 
1997) and be low in the optical; weak AGNs at last would break the spectral 
permanence by their non-equilibrium condition. This is supported by various 
lines of evidence: optical activity at very low levels is detected in a sizeable 
fraction of normal galaxies (Ho, Filippenko & Sargent 1997); X-ray galaxies with 
narrow optical lines are being detected, weak but so numerous as to conceivably 
saturate the XRB at ~ 10 kcV (sec Hasinger 1997); these advectivc accretion 
flows with low rj may constitute a stealthy addition to the relict BH masses. 

4. Population kinetics 

Bimodal fueling of BHs is conveniently described in terms of two components of 
the LFs at any L, z: 

N{L,z)=Ni{L,z)+N2{L,z) . (1) 

The component A'^i represents new BHs growing and flaring up during the host 
buildup, and dominates for M < 5 10^^ Mq or for z ^ 3 on average. Instead, 
represents BHs reactivated by interactions, and dominates in structures with 
M > 5 10^2 for 

z ^ 3, and especially in groups. 
The computation of A^i may be visualized in terms of the kinetic equation 
proposed by Cavaliere, Colafrancesco &: Scaramella 1991, which (setting the free 
/(M) from the integration so as to conserve the total condensed mass) generates 
the mass fuction N[M, z) of Press &l Schechter 1974: 

dtN = N/t+- N/t- with T- = 3t/2, r+ = r_(M/Mc)-("+2)/3 . (2) 

The positive driving term describes the build up of new host halos on a time 
scale tg ~ t. Their destruction on a similar time scale is described by the 
negative term. BHs masses and luminosities grow in the self-limited regime for 
At ^ rjtE ~ t mainly by such merging events rather than by continuos accretion 
onto a standing BH. Then the LFs are obtained in the form Ni{L, z) dL oc 
{At/t)pl{z)N{M,z)dM. 

A more complex story concerns the component N2 (see fig. 2). This is driven 
on the time scale r,. by the random reactivations of the A^^ dormant BHs, so the 
driving term for the kinetics is here A^/r^. Now the activity is supply-limited 
and mostly sub-Eddington (Cavaliere et al. 1988, Small 8z Blandford 1992). The 
BHs restart their bright career from a low luminosity Lj, and brighten up on the 
scale of the flyby time 2rg/V ^ fg/vg ~ 2 10^ yr, which is taken care of by the 
transport term 9i(LA2) in the kinetics. They attain an average L5, which sets 
the break in the LFs; statistically they quench off with increasing probability 

oc {LjL}^'^ Tgjvg (with ~ 1), and this sets the slope at the bright end. All 
that is expressed by the kinetic equation 
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Figure 2. An outline of the individual luminosities from intermittent 
interactions. Solid line: the average Lj, for accretion of gas in the host; 
dashed line: accretion of gas from satellites. 

To close the argument, the last term integrated over L provides one input 
to the number Nr of dormant BHs, the other coming from those still forming 
and flashing up as described by A^i; the component A''2 arises in groups from 
A^i and requires no independent normalization. We normalize A'^i to the data 
at 2 ~ 4. Adopting the stability limit Mbh oc {M/W^^MQf/''M (specifically, 
we use the coefficient lO'^ [(1 + z)/5]'^-^) has interesting virtues (HNR 1997): 
adequately flat LFs obtain at high-z; the number of QSs visible for At ^ 0.1 tE 
is consistent with that of the high-z, star-forming galaxies (Steidel et al. 1996). 

5. High and low-2; luminosity functions 

Fig. 3a shows the LFs for z ^ 3 computed in the critical universe using for the 
perturbations the tilted CDM spectrum normalized to COBE/DMR (see Bunn 
& White 1997); these are compared with optical data adopting the bolometric 
factor K = 10. For such redshifts not enough groups have yet formed for N2 to 
emerge to relevance; the evolution of Ni looks like the negative DE type. 

For z ^ 3, instead, N2 emerges and dominates the evolution. But to actually 
compute it, we have to consider where the main gas reservoir resides. Fig. 3b 
represents the low-z behavior when the accreted gas is provided mainly by the 
host reservoir (say, with a constant fraction used up in each interaction, also 
constituting the main gas sink); then —Mgas/Mgas — ~ —Lb/Li, holds. But 
in groups or clusters evolving hierarchically in the critical universe Tr oct closely 
obtains, since V (x t obtains for n ~ —2, while Ug oc Pu{z) oc t~'^ applies. The 
result is Li,{t) oc t^^°/'^^°; with r^o around 6 Gyr, scaled to groups from the classic 
census of local interacting galaxies by Toomre 1977, this reads oc (1 + z)'^. 

This implies LE dominant for 2; < 3, with LFs flattened at the faint end 
by the brightening over the flyby time, N{L) — > L-'^-'<'l^r(t) _ g^^. 
shows that quite some DE also occurs, due to the decreasing duty cycle of 
the reactivations governed by Tr{t) oc t. The LFs are steep at the bright end, 
N{L) exp [{— L / L),)'^ / <p]; but at low z the flat and slow Ni{L,z) remains 
exposed and flattens the overall slope, not unlike the data referenced in Sect. 1. 
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Figure 3. Optical luminosity functions {k = 10, = 10^^ erg s~^). 
Panel a from bottom: z = 6, 5, 4.3, 3.5. Panel 6 from top: z = 2.5, 
1, 0.5. Data: Kennefick, Djorgovski &; de Carvalho 1995; Schmidt, 
Schneider & Gunn 1995; Boyle, Shanks & Peterson 1988. 



Fig. 4b outlines the mass distribution of the relict BHs expected in most 
normal galaxies close to us, computed from Mbh = / dtL(t)/r)c^ with the lu- 
minosities evolving as above. 

When the gas in the host is depleted, the remaining reservoir is in satellite 
galaxies cannibalized out of an initial retinue gradually used up. With this 
prevailing, ~ const obtains with no LE (see fig. 2) , while the DE is enhanced 
as shown in fig. 4a. The gas available per event Mgas oc Mgat follows the dwarf 
galaxies distribution, to yield steeper LFs at the faint end. Wc expect such to 
be the case when M ^ 10~^Le/c^ holds, with ADAF prevailing; we relate this 
regime to the X-ray LFs of the NLXGs as given by Hasinger 1997. 



6. Conclusions and discussion 

We conclude that the coherent, non-monotonic QS evolution, with its scales of 
a few Gyrs and overall duration for ~ 12 Gyr, calls for a coordinating role of the 
surrounding structures to govern the accretion onto massive BHs. 

Remarkably, this is provided by the monotonic hierarchical cosmogony. 
Young host galaxies are assembled from subgalactic units, and then the ma- 
ture hosts are packed into larger and larger groups where they still evolve for 
a while by interactions. Both the assemblage and the interactions perturb the 
host gravitational potential and cause, besides starbursts, first rapid growth of 
BHs, then accretion episodes inevitably petering out in rate and strength. 

Thus the QS-galaxy connection is actually twofold. At low z much evidence 
relates AGN and QS hosts to interacting or to group environments. At high z 
one expects QSs associated with subgalactic star-forming blocks just looming 
out; one such instance may be the region singled out by Fontana et al. 1997. 



6 



-3 



-2 -1 



7 8 9 10 11 

M/Mg 



Figure 4. a) X-ray LFs from accretion of gas in satellite galaxies, 
with 1]/ K = 10~^; from top z = 1, 0.1. b) Mass distribution of the relict 
BHs, compared with data: masses from Kormendy &: Richstone 1995, 
and numbers from the bulge distribution after Franceschini et al. 1998. 

In this view the QS evolution should be marked by number increase as 
Mc{t) marches through the galactic range up to 5 10^^ M©; here the basic scale 
is tg. The turning point occurs as Mc{z) outgrows 5 IQ^"^ Mq at z ~ 3 ± 0.5; 
then the number, and even more the luminosities, decrease on the stretching 
scale Tr (X t. The LFs so computed are found to agree in some detail with the 
observations out to z ~ 5, specifically for the hierarchy provided by the tilted 
CDM perturbation spectrum in the critical universe. 

The predictions for z > 5 look rather dim; not many bright QSs are expected 
on the basis of the model successful at lower z (see fig. 3a). Actually these 
predictions sensitively depend on the two obvious parameters: the threshold for 
DM condensations in the Press &: Schechter mass function (in fig. 3a we use 
the canonical Sc = 1.69); the baryonic fraction Mbh/M packed in BHs (this 
decreases with M like M^/^, but is partially balanced by the factor (1 + z)^'^). 
The predictions for low z and weak L, instead, depend on whether the main gas 
reservoir is in the host or in satellite galaxies, as shown by figs. 3b and 4a. 

The quasar light history (of gravitational origin) computed from the above 
LFs is shown in fig. 1, and compared with the star light history (of thermonu- 
clear origin) given by Madau 1997; the latter at z ^ 1 is mainly contributed by 
the faint blue galaxies, which Cavaliere & Menci 1997 interpret as starbursts 
in dwarfs interacting in large-scale structures. Then the similarity of the two 
histories at such z reflects the basically similar run of the time scale oc 
for interactions in dense environments, in spite of two differences: the environ- 
ments arc constituted by condensing LSS or by virializcd groups, respectively; 
the FBG starbursts last longer than the nuclear activities. 

Looking briefly at other hierarchical cosmogonies, the hot -|- cold DM per- 
turbations - even when given all advantages like only 20% hot matter and suit- 
able amplitude to fit the data at z ~ 4 - produce at low z far too many bright 
QSs, reflecting the later collapses of galaxies in this version of the hierarchy. In 
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very open universes with standard CDM perturbations, we find for z < 2 too 
few bright QSs, due to group formation frozen long before and to gas reservoirs 
already exhausted. Low-density but flat universes fare much better, nearly as 
the critical, tilted CDM case. 
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